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INTRODUCTION  _ 

This  is  the  final  project  report.  fqr  AFOSR  grant  90-0307,  Wavelets 
and  Scattering.  The  original  goals  of  this  proposal  included: 

•  Obtaining  a  better  understanding  of  wavelets,  including  their  strengths 
and  weaknesses  in  areas  which  lie  in  the  mission  of  AFOSR. 

•  Applying  wavelets  to  electromagnetic  scattering  problems,  ultrasonic 
nondestructive  evaluation,  quantum  optics,  signal  analysis  and  the 
dispersive  structure  of  the  dielectric  response  of  water. 

Wavelets  are  very  promising  mathematical  and  computational  objects.  A 
single  function,  together  with  its  dilations  (translations  and  scale  changes), 
generates  a  basis  for  function  spaces  used  in  many  common  physical  situa¬ 
tions.  Scale  changes  provide  a  multiresolution  structure  which  clearly  and 
concisely  describes  elements  of  a  solution  or  operator  at  different  frequencies 
and  thereby  'unifies  diverse  areas  and  applications.  This  facilitates  the  trans¬ 
fer  of  methods  and  techniques  among  problems,  which  explains  the  appeal  of 
wavelet  methods.  Some  of  their  useful  properties  include  data  compression, 
noise  reduction,  edge  detection,  and  feature  extraction. 
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In  the  next  section.  Publications,  the  1994  publications  which  appeared, 
are  in  press,  or  submitted  are  listed  and  comments  are  given  on  what  has 
been  learned.  In  the  section  on  Follow  •  On  Projects,  the  next  studies  which 
are  made  possible  are  presented.  The  last  section  gives  the  Conclusion. 


PUBLICATIONS  1994 
Welland: 

Wl  G.  V.  Welland  and  M.  Lundberg  Construction  of  Compact  p- 
Wavelets,  Constructive  Approximation,  9,  347-370(1993). 

W2  R-  Torres  and  G.  V.  Welland,  The  Helmholtz  Equation  and  Trans¬ 
mission  Problems  with  Lipschitz  Interfaces,  Indiana  Univ. 
Math.  J.  42,  4,  1457-1485. 

W3  M.  Mitrea.  R.  Torres  and  G.  V.  Welland,  Regularity  and  approxi¬ 
mation  results  for  the  Maxwell  problem  on  Cl  and  Lipschitz 
domains,.  To  appear  in  the  proc.  of  conf.  on  Clifford  Algebras  held 
in  April  1993,  (12  typed  pp). 

W4  M.  Mitrea.  Torres  and  G.  V.  Welland,  Layer  potential  techniques 
in  electromagnetism,  submitted  to  J.  of  Int.  Eq.  and  Applic.  (29 

pp) 

W5  S.  K.  Bhatia,  V.  Lafcshminarayanan,  A.  Samal,  G.  V.  Welland  Pa¬ 
rameters  for  Human  Face  Recognition  submitted  to  J.  of  Visual 
Communicaton  and  Image  Representation 

(26  typed  pp) 

DeFado: 

Dl  C.  R.  Thompson  and  B.  DeFado,  Two-dimensional  image  analysis 
using  the  wavelet  transform,  in  Inverse  Problems  in  Scattering  and 
Imaging,  SPIE  1767,  Editor,  M.  Fiddy,  (SPIE,  Bellinghan  WA)  (1992) 
120  -  130. 

D2  B.  DeFado  and  S.-H.  Kim,  Non-uniqueness  in  direct  and  inverse 
electromagnetic  scattering  theory,  in  Inverse  Problems  in  Scat¬ 
tering  and  Imaging  SPIE  1767,  Editor,  M.  Fiddy,  (SPIE,  Bellingham 
WA,  1992),  pp  21-30. 
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D3  C.R.  Th^—^son  and  B.  DeFacio,  Information-to- noise  improve¬ 
ment  in  the  frequency  domain  using  the  wavelet  transform, 
in  Inverse  Problems  in  Scattering  and  Imaging,  1767,  Editor,  M.  Fiddy 
(SPEE,  Bellingham  WA,  1992)  pp  131  -  146. 

D4  D.M.  Patterson,  B.  DeFacio,  C.R.  Thompson  and  S.P.  Neal,  Wavelets 
and  their  applications  to  digital  signal  processing  in  ultra¬ 
sonic  NDE,  in  Rev.  Prog.  QNDE,  Edited  by  D.O.  Thompson  and 
D.E.  Chimenti  (Plenum,  New  York,  1993),  pp  719  -  726. 

D5  B.  DeFacio,  A.  Van  Nevel ,  and  0.  Brander,  Double  simple- harmonic- 
oscillator  formulation  of  the  thermal  equilibrium  of  a  fluid 
interacting  with  a  coherent  source  of  phonons,  International 
Workshop  on  Harmonic  Oscillators  NASA  Conf.  Proc.  1621,  Edited 
by  D.  Han,  Y.-S.  Kim  and  W.  Zachary  (NASA,  Greenbelt  MD,  1993) 
pp  309  -  322. 

D6  S.-H.  Kim,  G.  Vignale,  and  B.  DeFacio,  Frequency  and  wave- 
vector  dependent  dielectric  function  of  water-like  fluids,  Phys. 
Rev.  A  46,  7548-7560  (1992). 

D7  G.M.  D’Ariano  and  B.  DeFacio,  A  quantum  wavelet  for  quantum 
optics,  H  Nuovo  Cimento  B  108,  753-763  (1993). 

D8  S.-H.  Kim,  B.  DeFacio  and  G.  Vignale,  Refractive  index  of  water¬ 
like  fluids,  Phys.  Rev.  E.  48,  3172-3175  (1993). 

D9  B.  DeFacio,  Coherent-state  path-integrals  and  their  relation  to 
wavelets,  in  a  Festschrift  for  J.  R.  Klauder,  Edited  by  G.G.  Emch, 
G.  Hegerfeldt  and  L.  Streit  (Springer- Verlag,  New  York,  in  press)  18 

pp. 

DIO  S.-H.  Kim,  B.  DeFacio  and  G.  Vignale,  The  dynamic  dieleltric 
response  of  liquid  water,  submitted  to  Phys.  Rev.  E. 

Dll  B.  DeFacio,  S.-H.  Kim  and  A.  VanNevel,  Application  of  Squeezed 
States=Bogoliubov  transformations  to  the  statistical  mechan¬ 
ics  of  water  and  its  bubbles,  International  Workshop  on  Squeesed 
States  and  Uncertainty  Relations,  in  NASA  Conf.  Proc.  XXXX, 
edited  by  D.  Han,  and  M.  Rubin,  Y.  Shih  and  M.A.  ManTcov  (NASA, 
Greenbelt,  in  press)  13  pp. 
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D12  D.M.  Patterson  and  B.  DeFacio,  Wavelet  Inversion  of  Data  for 
Elastic  Wave  Nondestructive  Evaluation,  in  Inverse  Optics 
III,  SPIE  2241,  (SPIE,  Bellingham,  WA,  in  press),  13  pp. 

D13  H.  Kaiser,  K.  Hamcher,  R.  Knlasekere,  W.-T.  Lee,  J.F Anker,  B.  De¬ 
Facio,  P.  Miceli  and  D.  L.  Worchester,  Neutron  Optics  in  Layered 
Materials  in  Inverse  Optics  III,  SPIE  2241  (  SPIE  Bellingham 
WA,  in  press),  12  pp. 

FOLLOW-ON  PROJECTS 

•  Continuation  of  work  witb  (  McDonnell- Douglas  Corp.  CEM  group); 
Wavelet  methods,  matrix  sparsening  methods  for  3-D  elctromagnetic 
scattering.  (GVW) 

•  Statistical  mechanics  of  the  dialectric  response  of  water  (BDF  with  G. 
Vignale,  Cai  and  Welland). 

•  NDE  using  ultrasound,  especially  in  composites,  layered  materials, 
bond  adhesion  (with  Martin-Marietta,  Baltimore). 

•  p-biorthogonal  wavelet  project  with  B.  Jawerth  with  the  boundary 
adjustments  for  edge  sets.  (GVW) 

•  feature  extraction  using  wavelet  packets  for  the  development  of  a  face 
recognition  system  with  the  target  system  to  be  competitive  with  the 
human  capability  for  face  recognition. 

•  Novel  dielectrics  including  microwave  Penrose  tiling  structures,  wire 
structures,  periodic  arrays  of  spheres,  honeycomb  structures.  The 
band-gap  or  pseudo-gap,  switching  properties,  defect  and  the  disper¬ 
sion  will  be  studied  for  each  of  these  new  materials  (BDF  with  Satpa- 
thy). 

•  Inverse  Neutron  Optics  of  Layered  Materials  (BDF  with  H.  Kaiser  and 
D.  Worchester). 

CONCLUSION 

The  study  of  the  frequency  dependence  of  the  dielectric  response  of  water 
is  a  first  step  towards  a  quantitative  understanding  of  the  propagation  of 
electromagnetic  waves  in  tissue.  Such  understanding  is  required  to  answer 
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health  and  safety  questions  what  are  the  responsibility  of  the  Radiation 
Analysis  Division  of  the  Laboratory  at  Brooks  AFB  in  Texas.  The  squeezed 
states  developed  for  quantum  optics  will  be  studied  for  noise  reduction  in 
NDE.  It  is  expected  that  this  work  will  enhance  the  collaboration  of  BDF 
with  scientists  at  Martin-Marietta.  The  study  of  bubble  formation  and  their 
decay  addresses  questions  concerning  the  health  and  safety  of  ultrasound  in 
medicine,  since  these  decays  are  either  via  shock  waves,  or  initially  high 
velocity  jets  of  gas  and  matter. 

Numerical  work  based  on  multileveling,  sparsening,  and  conditioning 
methods  are  used  in  electromagnetic  scattering  problems  which  are  very 
difficult  to  resolve  because  of  technical  problems.  The  technical  problems 
derive  from  large  high  density  matrices  which  require  inversion.  Large  matri¬ 
ces  occur  in  cases  for  which  there  is  a  large  ratio  of  the  characteristic  length 
of  the  scattering  body  to  wavelength  of  an  illuminating  source  in  problems 
involving  body  element  methods.  This  is  at  the  heart  of  the  collaboration 
of  GVW  with  Lou  Mitschang,  D.-S.  Wang,  et  al.  at  McDonnell-Douglas 
(MDA).  This  has  been  extended  to  the  use  of  local  SVD  methods  and  the 
use  of  wavelet  packets  to  provide  a  controlled  sparsening.  The  goal  is  to  be 
able  to  solve  very  large  problems  using  sparsening  techniques  in  conjunction 
with  other  available  methods  or  methods  under  development,  while  main¬ 
taining  an  understanding  of  implications  with  respect  to  invertibility  and 
reliability  of  solutions  obtained  using  sparsening  techniques.  These  meth¬ 
ods  are  to  be  developed  with  the  objective  for  more  general  applicability. 
Other  groups  at  MDA  are  involved  in  this  project,  e.g.  the  CFD  group. 

Methods  using  wavelet  packets  are  being  developed  for  application  to  a 
machine  based  recognition  system  for  human  faces.  This  is  the  first  step 
in  a  project  to  obtain  good  representation  of  this  class  of  images.  One 
must  clearly  distinguish  recognition  from  identification.  The  goal  of  this 
effort  will  be  to  develop  a  system  that  rivals  the  power  of  the  human  system 
reported  in  the  paper.  Parameters  for  Human  Face  Recognition.  The 
methods  will  require  registration  of  images  through  shifts  and  rotations, 
which  requires  a  generalization  of  existing  methods.  This  problem  gathers 
together  all  of  the  ideas  in  the  present  project,  and  more. 

These  techniques  will  also  play  an  important  role  in  simulation  of  statis¬ 
tical  mechanics  models  and,  in  particular,  in  the  simulation  of  the  dielectric 
response  of  water  to  electromagnetic  pulses.  Our  studies  are  ongoing  and 
we  are  making  substantial  progress  toward  the  goals  stated  above. 

The  studies  completed,  so  far,  have  developed  methods  and  understand¬ 
ing  which  uniquely  qualify  the  Pi’s  for  these  follow-on  studies.  The  follow-on 
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studies  will  incorporate  refinement  of  analytic  modeling  of  physical  problems 
and  numerical  simulation. 

The  publications  and  preprints  listed  in  section  2  add  to  the  understand* 
ing  of  several  problems  of  interest  to  the  AFOSR.  The  p-wavelets  in  (Wl) 
provides  new  possibilities  for  segmenting,  compression  and  denoising  signals 
and  one  of  us  (GVW)  is  using  these  wavelets  to  study  edge  sets  with  Prof.  B. 
Jawerth.  The  papers  on  the  dispersion  of  water,  which  appear  in  references, 
[D6],  [D8],  [DIO]  and  [Dll]  show  established  structures  such  as  the  Debye  de¬ 
crease  in  dielectric  response  function,  and  such  potentially  new  effects  which 
have  not  yet  been  experimentally  observed  such  as  the  collective  dipolaron 
mode.  A  free  rotar  peak  was  found  in  the  far  infrared  frequency  region  at 
exactly  the  frequency  of  the  Simpson  et  al.  peak  which  was  experimentally 
reported,  but  had  not  been  obtained  in  other  theoretical  models.  However, 
the  peak  which  was  found  was  much  too  large  and  quantitative  agreement 
has  not  been  obtained.  The  dipolaron,  if  correct,  would  totally  change  the 
accepted  nature  of  the  relaxation  mechanism.  This  collected  mode  cannot 
be  ruled  out  by  the  data  which  exists  today. 

The  publication  listed  under  [D7]  presented  one  operator- valued  wavelet 
in  quantum  optics.  There  are  two  kinds  of  low-noise  quantum  states  in  op¬ 
tics,  coherent  states  and  squeezed  states,  and  it  is  well  known  no  additional 
such  objects  are  possible.  In  a  group  theoretic  approach,  the  coherent  state 
is  a  translation  in  Fock  space.  In  this  paper,  the  scale  change  of  wavelet 
theory  is  shown  to  be  implemented  by  the  squeeze  operator.  This  opens 
several  paths  for  study  of  noise  in  scattering,  in  addition  to  optics. 

The  publications  listed  in  [D4],  [D12]  and  [D13]  apply  wavelet  analysis 
to  inverse  scattering  of  elastic  waves  for  NDE  and  to  glancing  angle  ther¬ 
mal  neutrons  from  layered  materials  including  both  high  technology  metal 
composites  and  lipid  bilayers. 

We  take  the  opportunity  to  point  out  that  Brian  DeFado  was  elected  to 
Fellow  of  the  American  Physical  Society  in  November,  1993. 


We  have  not  applied  nor  obtained  any  patents.  We  have  made  no  in¬ 
ventions.  The  support  of  this  research  by  AFOSR  is  greatly  appreciated  by 
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Inverse  Problems  in  Scattering  and  Imaging 
SPIE  1767  -  4,  San  Diegj,  CA,  July  19-22,  1992 


Non-uniqueness  in  direct  and  inverse  electromagnetic  scattering  theory 


B.  De  Facio  and  S.-H.  Kim 

Missouri  University,  Department  of  Physics 
Columbia,  MO  65211 


Abstract 

General  statements  of  impossibility  can  be  important  in  science  and  en¬ 
gineering.  Ambiguities  in  inverse  problems  are  cases  of  non-uniqueness  where 
classes  of  different  objects  give  the  same  response.  A  strong  ambiguity  is  one 
which  no  additional  data  will  remove  the  non-uniqueness,  whereas  a  weak 
ambiguity  is  one  which  can  be  removed  by  additional  data.  In  direct  scat¬ 
tering  theory,  different  potentials  with  one  or  more  trapped  modes  may  give 
the  same  R(fc)  or  the  e“"R(fc)  where  a  is  a  real  parameter  at  all  wave-numbers 
Jfc.  In  three-dimensional  direct  scattering  theory,  different  material  media  and 
sources  J,  p  give  the  same  scattering  matrix  at  all  times  (or  wave  numbers) 
at  all  scattering  angles  and  all  incident  angles.  Examples  of  strong  ambigui¬ 
ties  will  be  given  including  one  where  a  temporal  relaxation  of  a  homogeneous 
body  is  equivalent  to  a  totally  different  time-independent  homogeneous  body. 
Weak  ambiguities  will  be  presented  including  both  examples  of  incident  scat¬ 
ters.  The  conditions  on  the  scatterers  at  spatial  infinity  and  their  trapped 
mode  bound-state  structure  will  be  given. 


1,  INTRODUCTION 


The  objective  of  inverse  problems1-4  is  to  study  families  of  models  for  solutions  in  & 
known  class,  which  are  amenable  computation.  Given  a  set  of  noise-free  data  the  question 
is  which,  if  any,  stable  solution  exists  in  the  class,  is  it  unique  and  how  it  can  be  calculated 
to  a  useful  accuracy  (in  a  reasonable  amount  of  time  on  an  available  computer).  If  the 
data  set  is  “too  small”  because  of  sampling  or  band-limiting,  or  if  the  signal  is  corrupted  by 
noise,  the  problem  becomes  ill-posed.  An  ill-posed  problem  has  a  solution  which  depends 
discontinously  on  the  initial  data,  so  that  it  is  also  unstable.  Most  of  the  interesting  applied 
inverse  problems  are  ill-posed.  At  the  1979  Delaware  Symposium  on  Ill-posed  Problems, 
P.C.  Sabatier  defined  interesting  inverse  problems  as  those  which  are  ill-posed. 

The  class  of  models  studied  in  this  paper  will  be  restricted  to  those  in  which  the  elec¬ 
tromagnetic  scattering  of  objects  is  described  by  Maxwell’s  equations  in  second  order  form. 
The  class  of  non-uniqueness  studied  here  was  first  discovered  by  Sabatier5-8  in  geophysical 
acoustic  scattering  and  generalized  to  electromagnetism  by  De  Facio*.  Corones  and  Winter10 
studied  the  related  homogenization  properties  of  a  scalar  electomagnetic  slab  (1-d)  problem 
using  splitting  methods.  These  ambiguities  clarify  the  proofs  of  uniqueness  for  electromag¬ 
netic  waves  which  are  found  in  refs.  (11  -  13).  Sabatier’s  work  was  focused  on  discontinuities 
at  interfaces  between  media.  Strong  dispersion  is  a  fundamental  property  of  the  dielectric 
response  of  biological  solutions14'15.  Inverse  problems  in  these  dispersive  media  were  reviewed 
in  ref.  (15).  The  physical  origins  of  dispersion  includes: 

(  1)  non-trivial  time  fluctuations  produce  frequency  u  dependence,  called  temporal  disper¬ 
sion,  upon  Fourier  transformation, 

(  2)  non-trivial  spatial  variation  in  an  object  gives  wave-number  k  dependence,  called 

spatial  dispersion,  upon  Fourier  transformation. 

It  is  the  temporal  dispersion  which  is  already  well  known  to  be  important  in  biomathematics 
and  at  present  the  spatial  dispersion  has  not  yet  been  established  experimentally.  It  will. 

The  identification  of  two  distinct  sets  of  coefficients  with  the  same  scattering  amplitude 
in  the  direct  problem  is  an  ambiguity  for  the  inverse  problem. 

Three  dimensional  inverse  problems  are  highly  overdetermined  and  the  exact  data  set 
which  is  to  be  used  requires  careful  consideration.  For  this  reason,  following  Sabatier3,  non¬ 
uniqueness  will  be  classified  as  a  strong  ambiguity  if  no  possible  additional  data  will  remove 
it  and  as  a  weak  ambiguity  if  there  is  possible  additional  information  which  will  remove  it. 
For  example,  near  field  information  of  the  scattered  field  would  remove  the  ambiguity  given 
in  ref.  (16)  but  not  that  of  ref.  (17).  If  it  is  possible  to  measure  the  near  field  response,  ref. 
(16)  is  a  weak  ambiguity  and  ref.  (17)  is  a  strong  ambiguity.  If  this  is  not  possible,  then  the 
objects  in  both  references  become  strong  ambiguities. 

It  is  also  possible  to  use  specific  examples  of  ambiguities  to  check  the  correctness  of 
electromagnetic  scattering  computer  programs,  as  pointed  out  by  Dr.  T.  Roberts18. 

The  organization  of  the  rest  of  this  paper  is  to  present  the  Theorectical  Framework  in 
Sec.  2,  Examples  of  the  homogenenization,  scaling  ambiguities  are  presented  in  Sec.  3  and 
the  Conclusions  are  given  in  Sec.  4.  The  Acknowledgements  and  References  follow  as  Sec. 
5  and  6,  respectively. 


2.  THEORETICAL  FRAMEWORK 


A  classical,  linear  electromagnetic  wave  at  frequencies  below  10u  Hz  satisfies  Maxwell’s 
equations  in  Minkowski  form1*: 


V  •£  =  />/,  (1®) 

7x£  =  -p  (16) 

V  •  J  =  0,  (lc) 

7x£  =  /e  +  p  (lrf) 


The  linear,  isotropic,  stationary  constitutive  relations  are: 

3  =  e£  =  f0£  +  P  (2a) 

5  +  fioM  (26) 

and 

If  =  o£  (2c) 


where  ( P ,  M)  are  the  electric  polarization  and  the  magnetization  in  matter  and  (e,  p,  a)  are 
the  permittivity,  permeability  and  conductivity  of  the  matter.  The  electromagnetic  field  pair 
(3,  3)  has  continuous  normal  componenets,  the  electric  field  £  has  continuous  tangential 
components  and  the  magnetic  intensity  £  has  coninuous  tangential  components  except  for 
the  idealized  (a  -»  oo)  a  perfect  conductor  where  the  discontinuity  in  £  at  the  interface  [£] 
has  the  non- zero  value 

(#]  =  [Jj]  =  £h  (3) 

where  £f  is  the  surface  current  density. 

The  model  to  be  studied  here*  is  an  electromagnetic  pulse  ( £0 ,  £o)  emitted  in  direction 
e  6  S2  in  a  linear,  isotropic,  homogeneous,  stationary  (LIHS)  medium  (l  toward  an  obstacle 
which  is  either  a  wavy  collection  of  layers  or  a  collection  of  compact  regions  fi2,  As,  •  •  ■ ,  fl/v+i 
where 

ft2  c  c  •••  c  n*+j  c  R3.  (4) 

The  constitutive  parameters  (€<,  /*„  a,)  and  shape  characterize  points  in  each  fl,.  Only  Qi 
is  LIHS  although  ft2  •  ■  ■Qw+1  will  be  taken  as  isotropic  and  stationary  for  simplicity.  The 
inhomogeneity  can  be  spatial  or  temporal  or  both.  If  the  electromagnetic  fields  and  sources 
(pi,  Jf)  are  in  certain  uniformly  continuous  Holder  spaces  given  in  Colton  and  Kress12  and 
used  in  ref.  (9),  then  the  scattering  amplitudes  all  exist.  The  scattered  fields  (£,  £)  are 
given  in  ref.  (9)  decomposed  into  pure  volume  integral  terms  were  called  diffuse  reflectors; 
surface  terms  integrated  over  discontinuities  in  [e],  [/i]  or  [o]  were  called  soft  reflectors;  and 
surface  terms  over  gradients  of  discontinuous  material  parameters  were  called  hard  reflectors. 
The  terminology  was  taken  from  geophysics.  One  consequence  of  this  decomposition  is  that 
a  fixed  frequency  pulse  cannot  distinguish  amoung  these  three  effects,  whereas  the  high 
frequency  part  of  a  wide- band  (narrow  time)  pulse  will  measure  only  the  hard  reflectors. 

It  is  well-known  that  variable-coefficients  partial  differential  equations  (pde’s),  without 
analytic  coefficients,  are  poorly  understood  even  in  their  existence,  uniqueness  and  stability 


properties.  The  simple  addition  of  spatial  or  temporal  dependence  adds  two  classes  of  pos¬ 
sible  non-uniqueness,  that  of  homogenization  and  that  of  Sabatier  scaling.  The  Sabatier 
scaling  may  be  discontinuously  or  dispersion  dependent  and  both  classes  can  give  either 
strong  or  weak  ambiguities. 

The  Sabatier  scaling  for  the  present  models*  involve  two  functions  (c„  p.)  which  can 
depend  on  t,  2  or  both.  However,  it  was  shown*  that  «,p,  =  1  is  required.  The  f-dependence 
is  particularly  important  for  water  and  biological  models  in  the  microwave  frequency  range 
and  has  been  the  subject  of  much  less  investigation  than  the  f-dependent  case.  For  the 
electromagnetic  fields  ( 3 ,  3,  3,  3)  and  their  sources  (p/,  //),  the  Sabatier  scaling  gives  a 
set  of  fields  (5~,  3~,  3~,  3~)  and  sources  (pj,  JJ)  defined  as 


&~  =  -D  =  ii.3, 

(5a) 

** 

II 

?  |- 

II 

•t: 

(56) 

&~  =  n.3, 

(5c) 

B-  =  P.5, 

(5  d) 

with  (,n,  —  1  (pointwise).  These  scaling  functions  correspond  to  the  multiplicative  inverses 
of  those  in  ref.  (9),  c,  (here)  =  1  /c,  (ref.  9),  etc.  The  constitutive  functions  (e,-,  p„  at)  and 
the  scaling  function  p.  can  depend  on  t  alone,  2  alone  or  t  and  2.  The  scaling  for  the  sources 
(p/,  Sj)  is  given  by 


and 


pH  =  P‘Pf  =  7^’ 

'I 

JJ  =  n.Jf. 


(ба) 

(бб) 


We  will  deviate  from  ref.  (9)  by  using  a  single  scaling  function,  which  one  of  us  (S-HK) 
found  to  exhibit  the  structure  in  a  clearer  manner.  Note  that  for  non-trivial  scaling  p,(.)  /  1 
in  the  interior  of  ft,  (i  /  1)  and  p,  =  1  (pointwise)  everywhere  in  fti . 

The  homogenization  occurs  when  (e„  m,  <r.)  are  functions  of  t  or  2  and 
post-homogenization  ambiguities  occur  when  the  scale  transformations  of  eqs.  (5a  -  d)  are 
applied  to  Maxwell’s  equations.  There  are  two  general  properties  of  the  Sabatier  scaling  of 
Maxwell’s  electromagnetism  which  will  be  presented  next. 

Lemma  1.  The  pairs  (Jf,  3)  and  (fy,  3")  have  the  same  conductivity  in  this  model,  for  all 
scaling  parameters  p,. 

Proof  Using  (5b)  and  (6b)  in  eq  (2c) 


J/  =  <t3  = «« Sy  —  <r(.3~ 
jj  =  o3~. 


Lemma  2  If  the  scaled  fields  (3~,  3~)  inside  ft,  (i  ^  1)  could  be  measured,  all  scaling  ambi¬ 
guities  would  be  weak  ambiguities. 


Proof  The  Poynting  vector  for  real  fields  (£,  8)  and  {S~,  8~)  are  given  by 


and 

Therefore, 


for  ft,  ft  ±1. 


S=£x8, 

S~  =  £~x  8~. 

S~  =  £(£  xd)  =  ft\§ 
*S 


Remarks 

(1)  On  page  865  of  ref.  (13),  a  correct  observation  is  made  that  different  conducting 
dielectrics  can  have  the  same  pointwise  Maxwell  equation  solutions.  The  authors  seem 
unaware  of  ref.  (9)  which  includes  this  result  as  a  special  case.  An  elementary  example 
will  be  given  in  the  next  section.  Lemma  1  shows  why  this  can  occur  structurally 
because  one  conductivity  can  correspond  to  different  pairs  ( If,  8)  and  (If,  £~)  which 
are  solutions  to  Maxwell’s  equations. 

(  2)  The  Poynting  vector  §  is  the  power  per  unit  area  flowing  in  the  direction  of  §  at  a 
point  f.  The  condition  §  ^  §~  shows  that  a  different  power  flows  inside  (i  /  1)  even 
though  the  external  fields  are  identical.  However,  this  lemma  is  inapplicable  to  remore 
sensing,  astronomy  and  many  objects  of  interest,  such  as  inside  a  nuclear  reactor  because 
observations  inside  the  object  are  impossible. 

In  the  next  section,  the  ambiguities  will  be  given. 


3-  AMBIGUITIES 

First,  the  homogenization  structure  of  Maxwell’s  equations  will  be  given  for  two  cases: 
0*2,  <2)  depend  only  on  time  and  then  (ft2,  e2)  depend  only  on  space.  The  conductivity  a  is 
carried  in  Sj.  The  space  dependent  case  has  been  far  more  frequently  studied  by  authors. 

Assuming  (e2,  fi2,  a2)  depend  only  on  time  t  and  taking  the  curl  of  eq.  (lb),  it  is 
straightforward  to  use  eqs.  (lc),  (Id)  to  obtain: 

Ai?  —  flltli  —  v(v  •  £)  +  9t(ft2ff)  +  10*2*2  -  +  (0tO*2*2))^]  +  0*2*2  +  2/i2f2)^  (7a) 

after  adding  a  background  term  -t\fi\£  to  both  sides.  Notations  of  £,  e2  are  used  for  partial 
derivatives  with  respect  to  time  of  a  single  vector  or  scalar,  dt((t2ff)  is  used  for  products  and 
A  is  the  Laplacian  operator.  A  similar  calculation  for  8  gives 


A 8  -  fi\(\8  =  v(V  •  8)  -  V  X  If  +  (0*2e2  -  ftiei)8  -  (d,(f t2e2))8]  +  ( ft2i2  +  2 it2f2)8.  (76) 


Next,  assume  that  («*,  *2,  <r2)  depend  only  on  space  f,  it  is  easy  to  take  the  curl  of  eqs. 
(lb),  (Id)  and  to  rearrange  into 


and 


a£  -  n\(iS  =  •  S)  +  it  if/  +  ((pa«a  -  ni(i)£\  +  (vmj)  x  8 


A&  -p  =  v(V‘&)  -  V  x  //  +  (/*2<2  -  ftifi)S-  v«2  x  £. 


(&») 

(86) 


The  coupling  between  £  and  £  in  eqs.  (8a),  (8b)  cannot  be  removed  for  general  (p2,  e2). 
At  fixed  frequency,  the  terms 


v(v  •  £)  =  (9a) 

«o 

v(v  B)  =  V  [^(VMa)  • 

=  (96) 

I/*?  .! 

where  eq.  (lc)  was  used  in  eq.  (9b)  and  pr  is  the  total  charge  density  (bound  plus  free)  can 
mix  with  additional  terms  to  provide  ambiguities.  It  is  clear  from  eqs.  (7)  and  (8)  that  fixed 
incident  (i.e.  scattering)  contains  little  information  on  the  details  of  the  obstacle  Sl2. 

For  a  metal  or  a  conducting  dielectric  the  time  dependent  (e2,  it2,  <r2)  case  the  right  hand 
side  is  either  equal  to  or  contains  a  term 

[Jiff  5]  =  3,0*2//) 

=  0*2*7/  -f  lt2Jj)nj  (10a) 

Example  1.  A  simple  calculation  shows  that  the  pairs  0*20,  //o)  and  (/*2i,  //1)  are  ambiguous, 
and  the  first  pair  gives  a  homogenization  of  the  second  when  the  directions  of  Jf0,  Jj  1  coincide 
and 


P20  =  1 
J/o  =  e-'W 

li2l  =  «-»«*«/* 

Jfl  =  (106) 

Calculating 

dt(ft2iJf\)  =  fclJfl  +  ft2\J/l 

=  ((ZM  +  ^)e-3PI«l/2]ePI«l/2 

=  =  dt(it20J/o). 

Therefore,  the  pairs  (p2,  J/)  =  (1,  e-<5|,|/2)  and  (e_w|t|/2,  e^|,|/l)  are  ambiguous  and  the  first  is 
homogeneous  in  time. 


Rsm&ikai 

The  decay  is  of  Debye  type  which  is  observed  in  polar  liquids  (including  H*0  and 
biofluids).  It  is  only  products  of  first  time  derivatives  of  /i2i  and  //,  which  must  have  decay 
faster  than  o(l/t)  at  large  times  for  scattering  theory  to  be  well-defined.  The  corresponding 
magnetic  field  3  behaves  consistently  since  y  x  Jc  is  perpendicular  to  hj  and  y  •  3  is  set  by 
continuity  of  the  tangential  fields.  This  gives  a  specific  example  of  the  ambiguity  mentioned 
in  ref.  (13). 

We  have  been  unable  to  construct  similar  examples  for  perfect  dielectrics  and  in  simple 
cases  such  as  Debye’s  exponential  time  dependence  times  a  polynomial  or  a  Guassian  time 
decay  times  another  polynomial  for  and  3  gives  uniqueness  when  the  uniform  Holder 
continuity  at  the  boundary  is  taken  into  account.  Since  this  occurs  in  such  limited  class  of 
dielectrics  and  fields  it  will  not  be  presented  here. 

The  examples  homogenization  from  spatially  dependent  constitutive  functions  are  espe¬ 
cially  anisotropic  because  of  terms  such  as  (y/i2)  and  (y«2).  Limited  angle  scattering  is  even 
more  questionable  than  usual  for  these  profiles.  Pure  metals  and  perfect  dielectrics  with 
only  spatial  dependence  with  continuous  boundaries  satisfy  Muller’s  uniqueness  theorm11 
and  this  can  be  relaxed  to  Lipschitz  continuous  at  the  boundary  with  exponent  o,  0  <  a  <  1, 
using  results  from  Colton  and  Kress12.  Hence,  any  fixed  angle  of  incidence  or  limited  angle 
ambiguity  is  weak  and  easily  removed. 

Now,  the  scaling  structure  will  be  given  by  substituting  eqs.  (5a  -  d)  and  eqs.  (6a,  b) 
into  Maxwell’s  equations.  Again  the  curl  of  the  curl  equations  will  be  rearranged  using  the 
assumption  that  (/i2,  «2)  depend  only  on  time  or  space  separately.  For  time  dependence  only, 
the  electric  field  satisfies  the  wave  equation 

A£~  —  =  v(y  •  3~)  +  JJ)  4-  0*2f2  —  +  [£2t2  +  2/j2c2  —  ^2e2^-)^~ 

+  [(£2*2  +/i2<2)  -  (£2*2  +  2/i2e2)—  —  ^2e2{^-i  —  2(^i)aV|fJ~  (11°) 

l  l  A1*  JJ 

and  the  corresponding  magnetic  field  satisfies  the  wave  equation 

A 3~  —  ftiCi3~  —  y(y  •  3~)  —  y  x  JJ1  +  (/i2e2  —  fii(i)3~  +  |/i2€2  +  2 £2e2  —  2/*2€2 ~| 3~ 

+  [(£2*2  +  £2*2)  -  0*2^2  +2£2«2)—  -  /i2e2{—  -  2(— (lit) 
L  K  !*•!*•  )  J 

The  structure  of  Maxwell’s  equations  gives  a  somewhat  simpler,  more  symmetric,  pair 
of  equations  for  the  case  of  a  spatial  dependent  scaling  function,  n,{2).  The  substitution  of 
the  scaling  equations  in  eqs.  (5),  (6)  into  Maxwell’s  equations  gives  the  pde  for  the  electric 
field  as 

-  V(V-^'“)+^2/7+(/i2f2-^iei)^"~-Vxl— (W.)x^]+[v/»2+- (V/*«)]xi?~.  (12°) 
The  pde  which  the  magnetic  field  satisfies  is 

&3~—n\ti3~  =  y(y-^'w)-yx/7+(/i2€2-/i,€j)^-V><(— (W.)x^~]-[yf2+— (y/i,))x£l~.  (126) 


The  coupling  between  the  electric  and  magnetic  fields  in  eqs.  (12a,  b)  is  identical  to  those 
in  eqs.  (8a,  b). 

The  time-dependent  scaling  case  is  richer  than  homogenization  case. 

Example  2  To  eliminate  the  dissipation  in  the  electric  field  £  in  the  general  case  of  time 
dependent  scaling  with  («*, pa)  time  dependent  only,  eq.  (11a)  requires  that 


A  2*2  +  2/ljfj  —  2/42*2 —  =  0. 

p* 

Since  n2t2  #  0,  it  can  be  divided  through  the  previous  equation  to  give 

+  2—  -  2—  =  0. 

P2  *2  p» 

This  can  be  integrated  directly  to  yield, 

P«(0  =  C3\/tl2(t)(2(t), 

which  is  the  most  general  scaling  allowed  for  («2,p 2)  dependent  on  time  only. 
Similarly  eq.  (lib)  is  dissipation-free  for  £  if 


^ 2*2  +  2pt2^2  ""  2^2*2 -  =  0 

which  can  be  solved  the  same  way  to  give  the  most  general  scaling  for  dissipation-free  £  as 

fi,(t)  =  C4y/(2(t)fi2{t). 


It  is  also  possible  to  have  a  weak  ambiguity  where  additional  measurements  can  restore 
uniqueness. 

The  space  -  dependent  scaling  case  is  complicated  and  coupled. 

Example  3  The  condition  of  dissipation- free  fields  when  (/i„/i2,e2)  depend  on  space  only  can 
be  solved  in  general.  From  eq.  (12b) 

~(V*2)  +  ~(VP«)  —  0 
*2  P« 

if  €2<,  ^  0  then  this  equation  can  be  integrated  directly  in  the  direction  of  the  gradients  to 
obtain 

*2  (*)P.(2)  =  Cl 

for  each  ci  €  Rl.  Hence,  there  are  ambiguities  for  dissipation-free  propagation. 

■ 


The  existence  of  trapped  modes  for  an  obstacle,  at  least  in  two  space  dimensions,  can 
be  described  by  repeatedly  performing  Darboux  transformations  on  a  scattering  operator 


for  the  object  without  any  interior  eigenvalue.  Each  transformation  produces  a  factor  of 
modulus  one  in  frequency  domain, 

'i=eiai.  (13a) 


k  +  iA 


k  -  iXj 


so  that  N  applications  gives 


n(B^)  - 


e*("i  +•■* +  OA/). 


(13*) 


In  terms  of  the  initial  reflection  operator  Ro(k),  the  same  obstacle  with  JV-trapped  interior 
modes  becomes 

RN(k)  =  c*(o‘  +  ••+  °'^)/20(Jk)  (13c) 

where  A,  is  the  wave  number  of  the  jth  trapped  mode.  The  Fourier  transform  of  eq.  (13c) 
with  ck  =  u>  gives  the  time-fomain  reflection  operator.  Thus  an  ambiguity  occurs  if  only  the 
modulus  squared  of  the  reflection  is  measured.  It  is  a  weak  ambiguity  for  a  classical  wave 
because  if  the  wave  amplitude  is  measured  the  phase  can  be  determined. 


4.  CONCLUSIONS 

The  non-uniqueness  of  the  time-domain  electromagnetic  scattering  problem  was  dis¬ 
cussed.  Two  kinds  of  non-uniqueness  for  inverse  problems,  a  stong  ambiguity  which  cannot 
be  removed  by  additional  measurements  and  a  weak  ambiguity  which  can  be  removed  by 
appropriate  additional  data  were  discussed  and  a  few  examples  were  given. 

In  electromagnetic  inverse  scattering  theory  a  recent  paper  by  on  of  us,  BDF#,  was  sim¬ 
plified  and  discussed.  Working  with  one  scaling  function  led  to  nicer  equations  but  two  were 
required  to  show  that  only  one  of  them  is  independent.  Two  general  features  involving  the 
Poynting  vector  of  scaled  interior  solution  and  another  showing  the  invariance  of  the 
conductivity  under  the  scaling  transformations.  This  actually  explains  an  observation  of 
Kabanikhin  and  Lorenzi15.  The  difference  in  ambiguities  due  to  different  non-constant  coef¬ 
ficient  functions,  called  homogenization,  and  ambiguities  due  to  scaling,  called  the  Sabatier 
scaling  or  electromagnetic  impedance,  were  derived  and  examples  were  presented.  Finally, 
these  examples  have  the  potential  to  act  as  checks  on  the  correctness  of  computer  programs 
for  electromagnetic  scattering  . 

It  would  be  interesting  to  use  new  technological  tools  such  as  molecular  beam  epitaxy 
(MBE)  to  use  these  ideas  to  attempt  to  fabricate  new  non-dissipative  materials.  Biological 
entities  may  already  be  using  the  idea  in  Examples  2  and  3. 

Electromagnetic  scattering  with  variable  permitivity,  permeability  and  conductivity  is  a 
rich  enough  model  that  the  surface  is  barely  scratched. 
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